(N 



in 
o 



BUCMP/02-04 
UMSNH-Phys/02-7 



Two- loop self-dual Euler-Heisenberg Lagrangians (II): 
Imaginary part and Borel analysis 



CN ' Gerald V. Dunne 

o : 

Department of Physics 
University of Connecticut 
Storrs CT 06269, USA 



Ph 

on 

os ; 

Christian Schubert 

(N 

Instituto de Fisica y Matemdticas 
Universidad Michoacana de San Nicolas de Hidalgo 



O ■ Apdo. Postal 2- 

o 

O 



J£V CP. 58040, Morelia, Michoacdn, Mexico 

■ schubert@itzel.ifm.umich.mx 



Center for Mathematical Physics, Mathematics Department 
Boston University, Boston, MA 02215, USA 

Qh! California Institute for Physics and Astrophysics 

Jj2 \ 366 Cambridge Ave., Palo Alto, CA 94306, USA 

> 

X 

Abstract 



We analyze the structure of the imaginary part of the two-loop Euler-Heisenberg 
QED effective Lagrangian for a constant self-dual background. The novel feature of 
the two-loop result, compared to one-loop, is that the prefactor of each exponential 
(instanton) term in the imaginary part has itself an asymptotic expansion. We also 
perform a high-precision test of Borel summation techniques applied to the weak-field 
expansion, and find that the Borel dispersion relations reproduce the full prefactor 
of the leading imaginary contribution. 



1 Introduction: imaginary part of effective Lagrangians 



The one-loop Euler-Heisenberg effective Lagrangians, for spinor QED |, §, | and scalar QED 
||], in a constant background electromagnetic field, have the following well-known integral rep- 
resentations: 
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Here T denotes the (Euclidean) proper-time of the loop fermion or scalar, and a, b are related 
to the two invar iants of the Maxwell field by a 2 - b 2 = B 2 - E 2 , ab = E • B. These La grangians 
have been renormalized [||, g |5| on-shell by subtracting, under the integral, the terms of zeroeth 
and second order in a, b. 



The effective Lagrangians ( |1.1| ) and ( |1.2[ ) are real for a purely magnetic field, while in the 
presence of an electric field there is an imaginary (absorptive) part, indicating the process of 
electron-positron (resp. scalar-antiscalar) pair creation by the field. The vacuum persistence 
amplitude is e lC , so the probability of pair production is approximately 2Im£ f|]. For example, 
in the case of a purely electric field, E, the effective Lagrangians (LI) and ( |1.2| ) have imaginary 
part given by: 
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where /3 = eE/m 2 . These expressions are clearly non-perturbative in terms of the field and 
coupling. 

The physical interpretation of these imaginary parts ( |1.3| ) and ( |1.4[ ) of the effective Lagrangians 
is that the coefficient of the k-th exponential can be directly identified with the rate for the 
coherent production of k pairs by the field |§ gj. This physical meaning of the individual terms 
of these series follows most directly from the alternative representation due to Nikishov |5|, 
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Here n p is the mean number of pairs produced by the field in the state with given momentum 
p and spin projection r, the =F refers to the spinor/scalar cases respectively, and p± is the 
momentum transverse to the field. An expansion of the logarithm in n p and term- by-term 
integration leads back to Schwinger's formulas 



]1.3| ) and (1.4). Thus the leading term in this 
expansion can be interpreted as the mean number n p of pairs in the unit 4-volume VT, while the 
higher {k > 2) terms describe the coherent creation of k pairs. This also explains the physical 
origin of the alternating signs in the scalar case (|bj) . 
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The exponential suppression factors in (Q and are very small for the field strengths which 
are presently possible for a macroscopic field in the laboratory. While positron production has 
recently been observed in an experiment involving electrons traversing the focus of a terawatt 
laser || , in this experiment the electric field strength is such that it is on the border-line between 
the perturbative and nonperturbative regimes, and so the positron production can be explained 
by multiphoton Compton scattering Q, rather than in terms of nonperturbative pair creation 
0]. This phenomenon has been analyzed in detail in ||, along with the prospects of observing 
nonperturbative pair creation using an X-ray free-electron laser which produces higher fields 
than those obtained by an optical laser. We mention also a recent analysis || which considers 
nonequilibrium and back-reaction effects. 

From a computational perspective, the imaginary parts of the effective Lagrangian can be com- 
puted in several ways. At one-loop, the most direct way is to consider the analytic properties 
of the integral representations in (|Ll]) and Ql,2| ). However, this approach is quite diffcult at the 
two- loop level for an electric field background, because the corresponding integral representa- 
tions are double integrals, with much more complicated integrands (see discussion below). An 
alternative approach is to exploit the well-known relation between the large-order divergence of 
perturbation theory and non-perturbative physics (l(], 11, 12]. For an electric field background, 



the weak- field expansion of the effective Lagrangian is a non- alternating divergent series, which 
signals the presence of an exponentially small imaginary part of the effective Lagrangian. Knowl- 
edge of the rate of divergence of the perturbation series can be used to deduce the exponentially 
small imaginary part of the effective Lagrangian, and therefore the pair production rate. 

This correspondence, between the divergent perturbative weak-field expansion and the non- 
perturbative imaginary part, is well understood at the one-loop level. For example, the weak- 



field expansions of the one-loop effective Lagrangians (LI) and (L2), specialized to the constant 
electric field case, are: 
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where the coefficients involve the Bernoulli numbers £> 2n , which alternate in sign and diverge 
factorially fast in magnitude [13, 14 1: 

S 2 n = (-ir +1 2^C(2n) (1.8) 

Here £( n ) denotes the Riemann zeta function |L 
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which is exponentially close to 1 for large n. Also noting that [13] 
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we see that the weak-field expansion coefficients in ( |1.6[ ) and ( |1.7| ) are 
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There is a precise one-to-one correspondence (as will be reviewed in section 3) between the 
terms in these large-order behaviors of the weak-field expansion coefficients, and the terms in 
the non-perturbative expressions ( |1.3| ) and ( |1.4j ) for the imaginary parts of the one-loop effective 
Lagrangians. 

At the two-loop level, the situation is even more interesting. The first radiative corrections to 
the Euler-Heisenberg Lagrangians ( |1.1| ) and ( |1.2| ), describing the effect of an additional photon 
exchange in the loop, were first studied in the seventies by Ritus |[^, [l(], [l?]]. Using the exact 
spinor and scalar propagators in a constant field found by Fock flg[| and Schwinger ||, and a 
proper-time cutoff as the UV regulator, Ritus obtained the two-loop contribution £^ in terms of 



a certain two-parameter integral. This integral representation was then used in [19] for deriving, 
by an analysis of the analyticity properties of the integrand, a representation for the imaginary 
part Im£g pin , analogous to Schwinger 's one-loop formula (|L~3|). Adding together the one-loop 
and the two-loop contributions, the imaginary part reads, in the purely electric case for spinor 
QED, 
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where a = §- is the fine-structure constant. The coefficient functions Kk(/3) appearing here 



were not obtained explicitly by [19]. However, it was shown that they have small (5 expansions 
of the following form: 
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Note that, for k > 2, these expansions start with terms that are singular in the limit of vanishing 
field (3 — ► 0, which seems to be at variance with the fact that these coefficients have a direct 



physical meaning. In [19] a physically intuitive solution was offered to this dilemma. Its basic 
assumption is that, if one would take into account all contributions from even higher loop 
orders to the prefactor of the k-th exponential, then one would find them to exponentiate in the 
following way, 
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Thus, it should be possible to absorb their effect completely into a field-dependent shift of the 
electron mass. Using just the lowest order coefficients in the small - /3 expansion of Kf.(^), those 
given explicitly in ( p.. 14 ), this mass shift reads 



m^(k,E) = m + -akckV eE — -akeE / m (1-16) 
19 1 these contributions to the mass shift have a simple meaning in the coherent 



As shown in [ 

tunneling picture: The negative term can be interpreted as the total Coulomb energy of attrac- 
tion between opposite charges in a coherent group; the positive one, which is present only in the 
case k > 2, represents the energy of repulsion between like charges. 
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It is important to note that this interpretation of the mass shift requires the mass m on the right 
hand side of ( |1.16| ) to be the physical renormalized electron mass of the vacuum theory. Only 
in this case the expansion of Kk{(3) has the form indicated in eqs.( |1.14| ). It is an interesting 
side result of the Lebedev-Ritus analysis that the physical electron mass can be recognized from 
an inspection of the two-loop effective Lagrangian alone, without ever considering the one-loop 
electron mass operator. 

Clearly, it would be of interest to understand in greater detail this prefactor series {(3) . It is not 
even clear from the work of Lebedev and Ritus whether or not this is a convergent expansion. 
In an effort to learn more about this prefactor, the imaginary part of the two- loop effective 
Lagrangian for a constant electric field background was studied using Borel techniques [pl|| . 
However, since no closed-form is known for the two-loop weak-field expansion coefficients in the 
case of a background constant electric field, this analysis was necessarily numerical. The leading, 
and first sub-leading, growth rates of the expansion coefficients were deduced numerically from 
the first 15 coefficients in the spinor QED case, using a brute-force expansion of the integral 
representation (a double integral) of the effective Lagrangian. This was sufficient to indicate 
that there are power-law corrections to the leading factorial growth rate, in contrast to the 
one-loop results (1.11) and ( |1.12| ) which only have exponential corrections to the leading growth 



rate. Using Borel dispersion relations, this indicates the presence of a prefactor series, as in 
the Ledebev- Ritus result ( |1.13| ). In this way, a numerical estimate was obtained for the next 
term in the expansion of the prefactor Ki(f3) of the leading exponential term (k = 1). This also 
provided an independent confirmation (albeit numerical) that the expansion proceeds in powers 
of \ffi. However, it was prohibitively difficult to go to much higher order, and so it was only 
possible to deduce one more term in the prefactor expansion. 

However, the case of a constant electric field background is not the simplest case one can study 



in this context. As shown in part I [22], for a constant self-dual Euclidean background, the 
two- loop effective Lagrangians, for both spinor and scalar QED, simplify dramatically, to such 
a degree that there are simple closed- form expressions in terms of the digamma function 
Such a constant background field satisfies 

"' " 2 



F/J.U = F\j,v = \e^ aP F^. (1.17) 



This has the consequence that 

F 2 = - fl (1.18) 
In Minkowski space, the self-duality condition ( [1.17 ) requires either E or B to be complex. 



This does not imply that such backgrounds are devoid of physical meaning; rather, the effective 



action in such a background should be understood in terms of helicity projections [£4], 22], and 
contains information on the photon amplitudes with all equal helicities (see part I). 

As shown in part I, for such a self-dual background, the two-loop spinor QED and scalar QED 
renormalized effective Lagrangians are: 
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as well as the important function 
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with rp being the digamma function ip(x) = T'(x)/T(x) fl3| , Oj. 

In [22], the parameter k, defined in ( |1.21| ), was taken to be real. In Minkowski space, this 
corresponds to taking the magnetic field B to be real, and the electric field E to be imaginary 
- so we refer to this case as the self-dual 'magnetic' case. In this case, the two-loop effective 



Lagrangians (1.19) and (1.20) are real. However, we could also analytically continue k — > in, 
which in the Minkowski language corresponds to taking E to be real and B to be imaginary. 
Thus, we refer to this case as the self-dual 'electric' case. For this background, with k imaginary, 
the two-loop effective Lagrangians ( |1 . 19 ) and ( 1.20 ) acquire exponentially small imaginary parts. 
As we shall see, these imaginary parts share many of the properties of the imaginary parts in 
(0) and (Q), for the case of a constant electric field background. 



The obvious advantage of the self-dual case is that since there are simple closed-form expressions 



for the two-loop effective Lagrangians ( 1.19 ) and ( 1.20 ), it is possible to study the imaginary 
part in much greater detail. In this paper, we analyze these imaginary parts, and compare the 
direct approach with the Borel dispersion relation approach. We are able to find the entire 
prefactor series for the imaginary parts. Since the spinor QED and scalar QED cases are so 
similar for a self-dual background, to avoid repetition, we concentrate on just one of them - the 
scalar QED case. In section 2 we present the details of the weak- and strong-field expansions of 
the one- and two-loop effective Lagrangians for scalar QED in a constant self-dual background. 
In section 3 we apply Borel summation techniques to the weak-field expansion, and use Borel 
dispersion relations to compute the imaginary part of the effective Lagrangian when the field 
strength is analytically continued k —* in. Section 4 contains some concluding remarks, and 
an Appendix contains the details of the derivation of the large-order behavior of the two-loop 
weak-field expansion coefficients. 



2 Weak-field and strong-field limits (scalar QED) 



In this section we discuss the weak- and strong-field limits of the on-shell renormalized one- 
and two-loop effective actions for scalar QED in a constant Euclidean self-dual background. We 
begin by recalling ( 1.2l| ) that k = "Vp where / is defined in ( 1-lcf) , so that "weak- field" 



2e v/ 2 ' 

means large k, and " strong- field" means small k. We first review the one-loop case, and then 
turn to the two-loop case. 



2.1 One Loop 



In the self-dual ('SD') case, the integral representation (1.2) for the renormalized one-loop ef- 
fective Lagrangian becomes 



m 



1 



(4vr) 2 4k 2 



dt 



-2nt 



sinh^t) 



(2.1) 



6 



using the Taylor expansion 14 1 



The weak-field (large k) expansion of this proper-time integral representation can be derived 
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where the expansion coefficients are (for n > 2): 
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The leading term in the weak-field expansion (^^) is 
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To derive the strong-field (small k) expansion, it is convenient to express (2.1) in terms of special 
functions as 
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where the function H(/c) is defined as pBl 
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and e'( — 1) ~ —0.16542. Note that this function is simply related to the function 
defined in ( |1.22 ), which appears in the two- loop expressions ( gig ) and (|L20|) by H'(k) = £(«)■ 

Thus, using the Taylor expansion of lnT(x), 
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the strong-field (small k) expansion of the one-loop effective action fl2.1|) is 
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The leading strong-field behavior in ( [2.9) ) is 
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Figure [T] shows a plot of (k 2 times) the exact one-loop effective action (1.12), compared with 
(k 2 times) the leading weak- and strong-field behaviors. (We multiply by the common factor 
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Table 1: The first 15 one- and two- loop weak- field expansion coefficients, and cffl appearing 



m and ( 2.16 ), together with their numerical magnitudes. Note that in each case the 

expansion coefficients alternate in sign and grow factorially fast in magnitude. Furthermore, the 
leading growth rate at one- and two-loop is the same. 



8 




Figure 1: This plot shows the exact one-loop scalar QED effective Lagrangian ([O]) [solid curve] 
in comparison with the leading weak-field (large n) [long-dashed curve] and strong-field (small 
k) [short-dashed curve] contributions ( [2.5] ) and ( 2.1 0| ) , respectively. We have multiplied through 
each function by a common factor of (An) 2 k 2 / m 4 for ease of presentation. 



of k 2 for clarity of presentation over a reasonable range of k). From Figure [l], it is clear that 



the leading behaviors ( |2.5| ) and ( 2.10| ) correctly capture the extremes, but fail to connect well 
in the intermediate region. (See [26[ for a simple approximate interpolation method to connect 
the weak- and strong-field extremes). 

From (|2.4| ) and (|1.8|) , it is clear that the weak- field expansion ( |2.3| ) is a divergent expansion. The 
first 15 expansion coefficients ch^ are tabulated in Table ^T|, together with their magnitudes. 
Notice that the coefficients alternate in sign (here k is real), and their magnitude eventually 
grows very fast. In fact, the leading growth rate of the coefficients at large n is factorial: 
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Figure ^ shows a plot comparing the exact one-loop result ( |2.1[) with successive truncations of 
the weak-field series expansion in ( |2.3| ). (Note once again that these plots have been multiplied 
by a common factor of k 2 ). We see that the successive truncations form an envelope around the 
exact result. This is typical behavior for an alternating asymptotic expansion. 

To probe more deeply the asymptotic nature of the weak-field expansion we consider the optimal 
asymptotic approximation p7| , |28| given by truncating the series in such a way that the error 
is minimized. (More sophisticated estimation techniques can be applied to the series, but the 
simplest one is sufficient here) . To find the optimal truncation point of an asymptotic series 

f{x) = J2(-V n CnX n (2.12) 

n 

we evaluate the terms c n x n and find the minimum value, say cn x n °. Then the optimal asymp- 
totic approximation consists of keeping terms up to n = Nq — 1. In this way the Nq 1 term, 
which is a measure of the truncation error, is minimized. Clearly, the truncation point depends 
on the value of expansion parameter x. 
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Nn = h + 7TK eXDf-r^— ) 
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% error 
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Table 2: This table shows the magnitudes of the terms cI^/k 2 ™ in the one-loop weak-field 
expansion ( |2.3| ), for four different values of k. In each column, note that the terms decrease in 
magnitude and then increase in magnitude. The order of the minimum term is n m i n , and this is 



compared with the estimate iVo of the order of the minimum term given in (2.14). The optimal 
partial sum is the partial sum up to the (n m j n — l) th term. This is compared with the exact 
value and the fractional error is shown as a percentage. 
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Figure 2: This plot shows the exact one-loop scalar QED effective Lagrangian ([O]) [solid curve] 
in comparison with successive partial sums of the weak- field (large k) expansion ( |2.3| ). The 
partial sum is shown for one term [long-dashed-curve] , two terms [medium-dashed curve], three 
terms [short-dashed curve], and four terms [dot-dash curve]. We have multiplied through each 
function by a common factor of (4ir) 2 k 2 / m 4 for ease of presentation. 



For the specific case of the weak field expansion ( |2.3j ), we define the remainder as 
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To find the optimal truncation point iVo, we minimize \Rn(k)\ with respect to N. The magnitude 
of this remainder is minimized when ln(2-7TK) = tp{2N — 1) ~ ln(2iV — 1) — l/(2(2iV — 1)) + .... 
Thus, a simple estimate of the point at which the remainder is minimized is given by 

N « - + ttk exp f - — J (2.14) 



2 \4ttk, 

In Table 2 we compare this optimal asymptotic approximation with the exact answer, for various 
values of k. Note that for k < 1 (i.e., "strong field", where / > ^-), the optimal truncation 
involves keeping just the first term in the weak-field expansion. Keeping more terms in the series 
makes the series estimate worse. 



2.2 Two Loop 

At two-loop, for scalar QED in a constant Euclidean selfdual background, the on-shell renor- 
malized effective Lagrangian is given by ( 1.20| ). The weak-field (large k) expansion of this result 
is 

^ SD, M = -(^E4?i (2-15) 
where the two- loop expansion coefficients are (for n > 2): 

fo\ 1 f 2fi — 3 „ 3 I B^k Son—Ik 1 / x 

c » , = (w{^ B2 "- 2+ 2£if(^%} <2 - 16) 
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Figure 3: This plot shows the exact two-loop scalar QED effective Lagrangian ( |1.20| ) [solid 
curve] in comparison with the leading weak-field (large k) [long-dashed curve] and strong-field 
(small k) [short-dashed curve] contributions ( |2.18D and ( |2.21| ), respectively. We have multiplied 
through each function by a common factor of (47r) 2 K 2 /(a7rm 4 ) for ease of presentation. 



This weak- field expansion follows directly from ( |1.20 ) and ( 1.22 ), together with the asymptotic 
(large x) expansion of the digamma function fl3| , |Tf| : 



ip(x) ~ lnx - ^- - E 



2x ^-J 2kx 2k 

k=l 



(2.17) 



The leading term in the weak- field expansion ( [2.15 ) is 



r (2)(SD), \ 



m 



(4vr) 2 128vr 2 « 4 



(2.18) 



The strong-field (small k) expansion follows directly from ( [1.20 ) and ( |1.22j ), together with the 
small argument expansion of the digamma function |13|, WM: 



-. oo 

T + Et- 1 )^)^- 1 

X k=2 



(2.19) 



This leads to: 



m 4 a 1 



(4vr) 3 k 2 



( 7 + ln«) -1 + -« + -k 2 ( 7 + ln«) - 3 £(-l) n C(n)« 



ra+l 



n=2 



— Li LXJ / Q \ Q OO Tl — 2 

+t k - E 9 c(») + (» + i)c(« + 1) ) « n + « E E (-i) B c(« - 0C(Q« B (2 

6 n=2 KZ 1 Z n=4 (=2 



20) 



The leading strong-field behavior is 



,(2)(SD) 



ma. 1 



scal (4tt) 3 k 2 



In k 7 



(2.21) 
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Figure || shows a plot of (k 2 times) the exact expression (|1.20f) , in comparison with (k 2 times) 



these leading weak-field and strong-field behaviors, (|2.1§| ) and ( [2. 21 ), respectively. This is 
completely analogous to the one- loop case (see Fig.Q). The extreme behavior for large and 
small k is well described, but they do not connect well in the intermediate region. 



The two- loop weak-field expansion ( |2.15 ) is a divergent asymptotic expansion; it is simply based 
on the asymptotic expansion ([2.17D of the digamma function iP(k). Figure || shows a comparison 



of the exact two-loop expression ( |l.20D with successive truncations of the weak-field expansion 
in (|2.15| ) (Once again, a common extra factor of k 2 is included for clarity of plotting). The 
pattern is much the same as in the one-loop case shown in Fig. |2l and is characteristic of an 
asymptotic expansion. 
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Figure 4: This plot shows the exact two-loop scalar QED effective Lagrangian (|L20|) [solid curve 
in comparison with successive partial sums of the weak-field (large k) expansion ( |2.15| ), The 
partial sum is shown for one term [long-dashed-curve] , two terms [medium-dashed curve], three 
terms [short-dashed curve], and four terms [dot-dash curve]. We have multiplied through each 
function by a common factor of (4ir) 2 k 2 / (anm 4 ) for ease of presentation. 



The two-loop expansion coefficients ( |2.16| ) alternate in sign and grow factorially fast in magni- 
tude, as shown in Table 1. In fact, the precise leading growth rate is 



.(2) 



, (-1)" 
'(2tt) 2 ™ 



T(2n - 1) 



1 + 



v 



(2.22) 



i 1 ) 



which is exactly the same as the growth rate of the one- loop cofficients ch' in ( 2.1 1| ) . This 
equality between the one-loop and two-loop leading growth rates can even be seen in the first 



15 coefficients shown in Table 2.1. This equality will become significant in the following when 
we discuss the Borel summation properties of the weak-field expansion. It also means that the 
estimate of the optimal truncation point is exactly the same as in the one-loop case, so we do 
not repeat the argument. 



3 Borel analysis 

The weak-field expansions Q and (glBD at one- and two- loop, respectively, are divergent 
asymptotic expansions. In this section we apply the technique of Borel summation to obtain 
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approximate expressions for the corresponding one- and two-loop effective Lagrangians. This 
enables us to test the Borel method with great accuracy because we are in the unusual position 
of having the explicit closed-form expressions (|2.1[ ) an d Ql.20| ) for the effective Lagrangians 
themselves, as well as the explicit expressions fly) and (|2.161) for the expansion coefficients to 
all orders. We also consider the issue of analytic continuation in the field strength, and compare 
the Borel results with the exact results. 

Strictly speaking, it is, of course, unnecessary to perform a Borel resummation in these cases, 
since we know the exact closed- form answers. However, such knowledge of the exact answer is 
rare, and our aim in this section is to probe the power, as well as the limitations, of the Borel 
approach in a case where a precise comparison can be made. 

Note of course that divergent behavior is not a bad thing; it is in fact a generic behavior 
in perturbation theory, as is illustrated in numerous and diverse examples in both quantum 
field theory and quantum mechanics |2^, 10, 11, 12 1. It is well known that knowledge of the 



divergence rate of high orders of perturbation theory can be used to extract information about 
non-perturbative decay and tunneling rates, thereby providing a bridge between perturbative 
and non-perturbative physics. 

To begin, we review very briefly some basics of Borel summation p0| , 57l [12^ , SlJ]. Consider an 
asymptotic series expansion of some function f(g) 



oo 

£ 

n=0 



a n g 



(3.1) 



where g — > + is a small dimensionless perturbation expansion parameter. In an extremely broad 



range of physics applications [11| one finds that perturbation theory leads not to a convergent 
series but to a divergent series in which the expansion coefficients a n have leading large-order 
behaviour 



On ~(-i)yr0in+ I /) 



(n — > oo) 



(3.2) 



for some real constants p, p, > 0, and v. When p > 0, the perturbative expansion coefficients 
a n alternate in sign and their magnitude grows factorially, just as in the weak-field expansions 
studied here: see ( [2.11 ) and ( 2.22 ). Borel summation is a useful approach to this case of a 
divergent, but alternating series. Non- alternating series must be treated somewhat differently 
(see below). 

To motivate the Borel approach, consider the classic example : a n = (—l) n p n n\, and p > 0. The 
series ([O]) is clearly divergent for any value of the expansion parameter g. Write 



fig) ~ £(-in w r / dte 

1 f°° , ( 1 A 

~ — / at I exp 

pg Jo \l + tj 



t 

pg 



(3.3) 



where we have formally interchanged the order of summation and integration. The final integral, 
which is convergent for all g > 0, is defined to be the sum of the divergent series. To be more 
precise ]30|, 27 1, the formula ( |3.3| ) should be read backwards: for g — ► + , we can use Laplace's 
method to make an asymptotic expansion of the integral, and we obtain the asymptotic series 
in fl3.1| ) with expansion coefficients a n = (— l) n p n n\. This example captures the essence of the 
Borel method. 
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For a non-alternating series, such as a n = p n n\, we need /(—<?). The particular Borel integral 
( |3.3| ) is |$(], 27] an analytic function of g in the cut g plane: |arg(g)| < n. So a dispersion 
relation (using the discontinuity across the cut along the negative g axis) can be used to define 
the imaginary part of f(g) for negative values of the expansion parameter: 



7T 



Im [f(-g)] ~ — exp[ 



(3.4) 



The imaginary contribution ( [3.4D is non-perturbative (it clearly does not have an expansion in 
positive powers of g) and has important physical consequences. Note that (3.4) is consistent 
with a principal parts prescription for the pole that appears on the s > axis if we make the 
formal manipulations as in (|3.3j): 



y p n n]g n^l_ f°° dt f M 

to 99 k Kl-tJ 



cxp 



t 

P9 



(3.5) 



Similar formal arguments can be applied to the case when the expansion coefficients have leading 
behaviour ( |3.2|) . Then the leading Borel approximation is 



f(9) 



1 



dt 



/x Jo t \l + t 



1 



v/fi 



exp 



\pgJ 



(3.6) 



For the corresponding non-alternating case, when 5 is negative, the leading imaginary contribu- 
tion is 



Im [/(-<?)] 



7T 



/ 1 y/n 


1 1 \ 


— J exp 




\pgJ 


\P0/ 



(3.7) 



Note the separate meanings of the parameters p, /i and that appear in the formula fl3.2j ) for 
the leading large-order growth of the expansion coefficients. The constant p clearly combines 
with g as an effective expansion parameter. The power of the exponent in fl3,6j ) and ([O]) is 
determined by fi, while the power of the prefactor in (3.6) (|3.7| ) is determined by the ratio — . 



To illustrate this correspondence, between the perturbative expansion coefficients (|3^) and the 



imaginary part (3.7), it is now a straightforward exercise to check that for the one-loop electric 
field background, the weak- field expansion coefficients in ( |1.11| ) and ( 1.12[) correspond to the 
imaginary part of the respective effective Lagrangians given in ( |1.3[) and (1.4). 

These formulas ( |3.6j ) and ( |3.7| ) are somewhat formal, as they are based on assumed analyticity 
properties of the function f(g). These can be rigorously proved when the expansion coefficients 
are precisely given by the gamma function form in ( |3.2[ ). But things are less clear when this 
expression gives just the approximate large order behavior of the coefficients. The Borel disper- 
sion relations could be complicated by the appearance of additional poles and/or cuts in the 
complex g plane, signalling new physics, such as, for example, renormalons |3l], |3^]. Also, the 
expression ([T7]) assumes a principal parts prescription for the poles. Our analysis below will 
show that these assumptions are justified for the one- and two- loop effective Lagrangians (|2.l| ) 
and (|L20|) for scalar QED in a constant Euclidean self-dual background. The same is true for 
the spinor case, but we do not repeat the analysis here: it is very similar. 
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3.1 Borel Analysis of One Loop Case (Scalar QED) 



To begin, we write the one-loop weak-field expansion coefficients Cn^ in ( |2.4|) in a form which 
we can easily relate to This is clearly already true of the leading behavior ( |2,ll|) , but 

we can also quantify the corrections. Recall Euler's formula ( |1.8| ) that relates the Bernoulli 
numbers to the Riemann zeta function. We shift the index of the expansion coefficients so that 
the weak-field summation in (2J3) begins at n = 0, as in (3T). Thus, we define new expansion 
coefficients 



2n+4 



(2n + 2)(2n + 4) 
(-l) n 



(27T) 



2ra+4 



r(2n + 3) +T(2n + 2)] 1 + 



1 



+ 



1 



22n+4 1 32n+4 



+ 



(3.8) 



Keeping the exponentially leading term 1 in the curly bracketed factor, we see that the coefficient 
On is the sum of two terms, each of which has the form in (|3.2j), with jjL = 2, p = l/(27r) 2 , and 
v = 3 and 2 respectively. So, a direct application of the Borel integral expression ([O]) leads to 
the leading Borel approximation 



C 



(i) 

leading 



(«) 



m 



1 



(4vr) 2 2tt 4 k 2 



,-2/tt 



dt 



1 + 



-(t + 2k t 2 



(3.9) 



Figure |5| shows a comparison of this leading Borel approximation with the exact effective La- 
grangian (as before, each with an extra factor of k 2 ). Note that it does very well at large n 
(weak field), but less well at small k (strong field). Nevertheless, it does much better in the 
strong field regime than the partial sums of the weak-field expansion (compare with Fig. |2|). 
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Figure 5: This plot shows the exact one-loop scalar QED effective Lagrangian (^^) [solid curve] 
in comparison with the leading Borel approximation in (]3.9|) [dashed curve] . We have multiplied 
through each function by a common factor of (4-7r) 2 K 2 /m 4 for ease of presentation. Note that 
the agreement is much better than for the weak-field (large k) expansion shown in Figs. 1 and 
2. 

Successive subleading corrections to this leading Borel approximation follow from considering 
the (exponential) corrections inside the curly brackets in the expression (|3.8| ) for the expansion 
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coefficients. Once again, these terms are of the "standard form" in (3.2), now with p = l/(2Tvk) 



for k = 2, 3, . . .. Thus, each term can be resummed using the Borel formula (3.6). This leads to 



4Ii(«) = 0^2 Eu * T^H' + 2 ^ 2 ) ( 3 - 10 ) 



m 4 -I 00 -i roo p -2nt 
(4vr) 2 2vr 4 K 2 2_, fc 4 J Q 1 + 

Keeping successive terms in the expansion (fUOp rapidly improves the agreement beyond that 
shown in Figure |B| for the leading Borel approximation ( |3.9| ) . 

In fact, the expression ( |3.10| ) could have been derived directly by interchanging the summation 
and integration. To see this, note that by an integration by parts we can rewrite the proper-time 



form (2A) of the effective Lagrangian as 



In the large k limit we can use the small t expansion of the coth function [13| 



coth t = tV 3.12) 

from which the Borel result fl3.10| ) follows. Thus, the proper-time integral representation fl2.1| ) 
is precisely the Borel integral representation of the divergent asymptotic weak-field series ( |2.3| ). 
And, conversely, the weak-field series expansion fl2,3|) is just the large k asymptotic expansion 
of the proper-time integral expression (p.l|). This is completely analogous to what happens in 
the case of a constant magnetic field background 

Now we turn to the Borel analysis of the imaginary part of the effective Lagrangian. Recall 
that in the magnetic/electric case, the transition from a magnetic to an electric background 
involved the perturbative replacement B 2 — > — E 2 , or B — ► iE. This is because the Lorentz 
invariant combination is (B 2 — E 2 ). Similarly, in the self-dual case, we change from the self-dual 
'magnetic' case to the self-dual 'electric' case by the the analytic continuation k — > ik. Under 
this transformation, the only change in the weak-field expansion ( [2 .3D is that it now becomes a 
non- alternating divergent series. Thus, according to the Borel dispersion relation formula ( |3.7| ), 
the effective Lagrangian acquires an imaginary part. Using this Borel dispersion relation formula 
( |3.7| ), together with the explicit expression ( |3.8[) for the expansion coefficients (but now with the 
alternating (— l) n factor suppressed), we find 



Im 



This agrees precisely with what one finds by a direct analytic continuation of the proper-time 



integral representation (2.1), together with a principal parts prescription for the poles. 



3.2 Borel Analysis of Two-Loop Case (Scalar QED) 



The Borel analysis at two-loop is more complicated because the corrections to the leading be- 
havior ( 2.22 ) of the two-loop weak-field expansion coefficients ( 2.16| ) are considerably more 
complicated than for the one-loop case ( |3.8[ ) where the corrections were simple exponentials. 
Nevertheless, an analysis of the large order behavior of the two-loop coefficients ( |2.16D uncovers 
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a remarkably rich number theoretic structure, and permits a Borel analysis of the complete 
prefactor to the leading exponential behavior. This will be shown in this section. 

As before, we shift the summation index in the weak-field expansion (EUHD so that it be gins at 
n = 0. So we define the expansion coefficients 



(2) = J2) 



-n+2 



(3.14) 



where ch^ are the two-loop weak-field expansion coefficients given by (2.16). In the Appendix 
we show that if we neglect sub-leading exponential corrections (this is analogous to just taking 
the 1 term in the curly brackets in ( |3.8| ) in the one-loop case), thenQ as n — > oo, 



,( 2 ) 



-l) n 



(27T) 



2n+4 



r(2n + 3)- i r(2n + 2)--— 



,|''-" l yV t 4-ai' 



Z=2 



2/ 



(3.15) 



This should be compared with the analogous formula for the one-loop case when we make the 
same approximation of ignoring the exponential corrections in (|3.8|): 



-l) n 



(27T) 



2ra+4 



[r(2n + 3) + T(2n + 2)] 



(3.16) 



In each case we have written the expansion coefficients as a sum of gamma functions, each term 
of which has the general form in ( |3.2|) . The difference is that at one-loop there are only two 
terms, while at two-loop there is an infinite series of terms. This translates directly into the fact 
that in the one-loop imaginary case, the prefactor of the exponential in ( 3. 13] ) has two terms, 
while in the two-loop case we will find an infinite series of terms in the prefactor (see ( |3.20| ) 
below). 

Given the expansion ( 3.15| ) it is straightforward to use the Borel formula ( |3.6| ) to write the 
leading Borel approximate integral representation for the two-loop effective Lagrangian: 



Aeading ( K ) 



m 



an 



(Air) 2 7T 4 K Jo 



dt 



e~ 2K H 2 
1+5" 



J_ _ 6Kt ^ {-l) l B 2 i 
8«t {2nf y 21(f)* 



1=2 



(3.17) 



Using the expansion |13J : 



Re[ip(l + iy)] =\ny + ^T 



(-iy- i B 2l 



i=i 



2ly 



21 



(3.18) 



we can rewrite the leading Borel approximation ( [3.17 ) as 
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leading 
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an 
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dt 



e -2nt t 2 
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2kI 



6nt 



Re v(i + — : 

7T 



In — 



(3.19) 



In Figure ^| we compare this leading Borel integral representation ( 3.19| ) with the exact two-loop 
effective Lagrangian ( |1.20D . Notice that the agreement is not as good as in the one-loop case 

The derivation of this expansion requires some nontrivial results in number theory. This is explained in the 
appendix. 
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Figure 6: This plot shows the exact two-loop scalar QED effective Lagrangian ( |1.20| ) [solid 
curve] in comparison with the leading Borel approximation in ( |3.19| ) [dashed curve]. We have 
multiplied through each function by a common factor of (4ir) 2 k 2 / m 4 for ease of presentation. 
Note that the agreement is not as good in the small k (i.e., strong field) region as that obtained 
in the one-loop case in Fig. 0. See text for an explanation. 



(compare with Figure |5|). This can be traced to the fact that the second term in the expansion 
( 3. 15] ) of the two- loop weak- field expansion coefficients is negative. On the other hand, the 
second term in the expansion ( 3.1G| ) of the one-loop weak-field expansion coefficients is positive. 
This second term, in each case, leads to a log divergence as k — > 0. In the one-loop case this 
log divergence has the same sign as the log divergence of the full answer (see ( |2.10|) ), but in 
the two-loop case the log divergence coming from the leading Borel expression (|3.19|) has the 
opposite sign compared to the log divergence of the exact answer in the strong field (small k) 
limit - compare with (|2.21| ). Hence, this mis-match becomes more dramatic at small 
clearly seen in Figure ||. 

Now consider the analytic continuation k — > in. The only difference in the weak-field expansion 
is that the series (2.15) becomes a non-alternating series. Thus, according to the general Borel 
dispersion relation, there should be a series of imaginary parts, which can be deduced from the 
formula (|3.7| ) . This leads immediately to (here we use a principal parts prescription for the pole 
at s = ir) 



Im 



^leading 



m 
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air 



(4vr) 2 2k 
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2ttk 2tt 



3k™(-1) 1 B 21 



i=i 



21k 



21 



-2ttk 



(3.20) 



Note that the dominant part, at large k, of the prefactor (the 1 in the square brackets) is 
exactly the same as the one-loop leading term in fl3,13|) , except for an overall factor of an. In 
the language of the Borel analysis, this can be traced directly to the fact that the leading growth 
rates of the one- loop and two-loop expansion coefficients in ( 2.11 ) and ( |2.22 ), respectively, agree 
precisely. (The overall factor of air in the two-loop case was explicitly separated out in front 
of the two-loop weak-field expansion (^15|)). This is the same behavior as is found in the case 
of a constant electric field background, where the leading two-loop imaginary part is air times 
the leading one-loop imaginary part [compare Equations (b4), ( 1.13 ) and ( Qj )], and which is 
reflected also in the same leading growth rates of the one- and two-loop weak-field expansion 



coefficients [21|. 
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We can now compare this with the imaginary part obtained from the exact closed-form expression 
(|L2C| ). Using the fact that [l3|] 



1 7T 

Im = h — coth(7TK) 

2k 2 



(3.21) 



we deduce that 
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where we have defined the shorthand function 
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(3.22) 



ip(x) = ip(x) — Inx + 



2a- 



(3.23) 



Thus, the leading Borel expression ( |3.20| ) is simply the k = 1 term (i.e., the leading exponential 
contribution) in this expansion. This is completely consistent with our earlier approximation 
of only keeping power-law (but not exponential) corrections to the expression for the expansion 



coefficients. So we conclude that the direct use of the Borel dispersion relation (3.7), together 
with the implied use of a principal parts prescription, is consistent in the two-loop case also. 

It is interesting to see that in ( |3.22 ) we have the complete Lebedev-Ritus prefactor expansion, for 
any instanton index k. This should be compared with the two-loop electric field case ( 1.13| ) and 
( 1.14| ), where only the first few terms of this prefactor are known. In particular, it is now clear 
that this expansion is itself a divergent asymptotic series, a fact which was left undecided by the 
Lebedev-Ritus analysis. Another interesting issue is the following: comparing the self-dual case 
result ( 3.22j ) with the corresponding one-loop result ( |3.13| ) we see that, for any non-zero 'field 
strength' /, the k-th prefactor at two-loops will dominate over the corresponding one-loop term 
for sufficiently large k. On the other hand, in the self-dual case, it is straightforward to verify that 
the total two-loop contribution ( 3.22j ), when one sums over all k, is always smaller (essentially 
by a factor of air) than the total one- loop contribution ( |3.13| ). It would be very interesting to 
see if a similar thing happens in the electric background case, where the individual k terms in 
the sum have a direct physical meaning in terms of a coherent multi-pair production process 0] . 



4 Conclusions 

To conclude, we have presented an analysis of the imaginary part of the one-loop and two-loop 
effective Lagangians for scalar QED in a constant self-dual background. At both one- and two- 
loop, these self-dual effective Lagrangians have many properties in common with the effective 
Lagrangians for a constant magnetic or electric field background. They have similar weak- and 
strong-field expansions. Also, the correspondence between one- and two-loops is very similar in 
the magnetic/electric and self-dual backgrounds. In particular, in the self-dual case we found 
that the imaginary part of the two-loop effective Lagrangian has a leading term equal to the 
leading term of the imaginary part of the one-loop effective Lagrangian, multiplied by an overall 
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factor of air. In terms of the weak- field expansion coefficients, this is reflected in the fact that 
the one-loop and two-loop expansion coefficients (with a factor of onr extracted in the two-loop 
case) have identical leading growth rates at large order of perturbation theory. We stress that 
this agreement between leading growth rates is highly nontrivial - it only occurs if one uses the 
consistently renormalized mass parameter in the expansion, and thus it is sensitive to the finite 
part of the mass renormalization that enters this analysis at two-loops. These facts also hold 
for the constant electric field background. It is also consistent with the exponentiation factor, 
exp(a7r), found by Affleck et al for the case of pair production in a weak electric field at strong 



coupling in scalar QED [34]. Indeed, based on the expectation that this exponentiation occurs 
also in the self-dual case, we conjecture that at loop order I, the weak-field expansion will take 
the form 

L 1 > (1-1)1 (47r)2 ^ 

where for each I, the expansion coefficients c4 have the same leading large n growth rate as oh . 



Thus, for scalar QED, the c4 should grow factorially exactly as in (2.11), while for spinor QED 



the leading large n growth should have an additional factor of —2. This also has implications 
for the leading large n behavior of the zero momentum limit of the all "+" helicity amplitudes, 



as is explained in [22] at one- and two-loop. In particular, this conjecture (4J) for the growth 
with loop order can be immediately extended to the low energy limit of the 'all +' N-photon 
amplitude at large N. At first sight this might seem a puzzling result, since it suggests that 
the N-photon amplitude, if only in one helicity component and in the low energy limit, can 
be expressed as a convergent power series in a for sufficiently large N. This runs counter to 
expectations based on well-known general arguments which indicate that the loop expansion 
should yield an aymptotic series. This apparent contradiction might be resolved either by the 
appearance of poles in the complex a plane at higher loop orders, invalidating the naive use 
of the Borel dispersion relations, or, what seems much more likely, by a slowing down of the 
convergence in N of the ratio f^mrj with increasing I. Such a non-uniformity of the convergence 
might be visible already at the three-loop level, which lends additional motivation to pursuing 
the calculations presented here to even higher loop order. 

We also believe that the results we obtain here for subleading behavior, in the self-dual back- 
ground, can be taken as illustrative of the constant magnetic/electric case also. The advantage, 
however, of the self-dual case is that there exist closed- form expressions for the two-loop renor- 
malized effective Lagrangians. This means that the analysis is much more complete in this case. 
We were able to find the full imaginary part in the self-dual 'electric' case, and we showed that 
this is consistent with what one obtains by a Borel dispersion relation treatment of the weak- 
field expansion, which is a non-alternating divergent series. At two-loop, the prefactor of each 
exponential term in the imaginary part, has itself an asymptotic expansion. 
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A Asymptotic expansion of two-loop coefficients in scalar QED 



(2) 

In this Appendix we compute the large n - expansion of the coefficients c n , appearing in the 
two- loop weak field expansion for scalar QED, up to subexponential terms. We recall from 
(|2.15| ) that the expansion of the final two-loop formula (|l.20|) yields the following expression for 
the coefficients: 
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An alternative form for these expansion coefficients arises when one writes ^-^Jf^ (k) as a single 
proper-time integral (see Equation (3.17) in part I p2[): 
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This second form is more convenient for studying the large order behaviour of c n ■ The equiv- 
alence of the formulas ( |A.1| ) and ( A.2 ) is not immediately obvious. However, it can be shown 
using Euler's identity for Bernoulli numbers 
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together with a more involved identity known as Miki's identity |37, 38]: 
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It is the second form, in (|A.2[) , for c n which we will use in the following for obtaining the 
large n behavior of these coefficients. To start with, we note that the large n expansion of the 



identities for sums of products of Bernoulli numbers appear in many contexts in physics and mathematics 
il|36|. 
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first term in the curly brackets is straightforward, using Euler's formula (1.8) which relates the 
Bernoulli numbers to the zeta function, which is of order 1, with exponentially small corrections. 
The difficulty is in the last term in ( A.2| ), because the index on each of the Bernoulli numbers 
need not be large just because n is large. Thus this last term, which involves the finite sum of 
a product of Bernoulli numbers must be modified somehow to bring it into a more manageable 
form. We therefore use the Euler identity (|A.3|) to rewrite (for n > 2) 



n-l 



2n-2fc 




n-l 




+2(-l)' 



(gn - 1)! 

(27T) 



B>2kE>2n-2k 

" _1 ((2k)((2n - 2k) 



2)1 



E 

k=l 



n 



k 



2n(2n - 1) 



+ 2(-l) n 



[2n 



(2tt) 



1)! g Q{2k)Q{2n - 2k) 



2n 



k=l 



k 



(A.5) 



The advantage of these manipulations is that we can now separate out the leading large n growth 
rate of the coefficients and write: 
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In the approximation which neglects exponentially suppressed terms, the large n behavior of 
the first two terms is simple, since the £ - factors can be replaced by unity. The large order 



behavior of the tp function follows from using ip(m) = ip(m — 1) + 
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m— 1 ' 



and the asymptotic 



expansion(2.17). Thus, we find (here and in the following, ~ denotes equality up to terms which 
are exponentially small for large n) 
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Note that it will turn out to be convenient to perform the expansion in n — 1 rather than in n. 



The third term in (A.6) is the problematic one. It involves a folded sum of zeta functions. While 
in many cases such sums can be reexpressed in a form which involves the zeta function only 



linearly 



], for the case at hand no such formula appears to be known. However, we can 



still obtain its asymptotic expansion from Euler's formula 
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We can immediately discard all terms where not either r = 1 or s = 1, since they are exponen- 
tially small. Thus we have to consider the following three contributions: 
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In the last step we used the combinatorial "normal ordering" identity 
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(m > 0), where Sm ' denotes the Stirling number of the second kind, which is the number of 
ways of partitioning a set of m elements into k non-empty subsets. Explicitly: 
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The s - sum can be done in closed form after using a partial fraction decomposition to write 
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Thus, again using ( p. 17 ) we obtain the asymptotic expansion of Z n : 
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where the are the numbers defined in ( |A.14 ). Finally, putting everything together, the 

(2) 

final result for the asymptotic expansion of the coefficients Cn becomes, up to sub exponential 
corrections, 
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Observe that the leading 0(log(n — 1)) and 0((n — 1)°) terms cancel out in the sum. This is 
a nontrivial self-consistency check. We note that, for the constant term, the occurrence of this 
cancelation depends on the value of the finite part of the mass shift in the mass renormalization 
[22]; it happens only if the renormalized scalar mass is the physical one. That the leading 
asymptotic behaviour depends qualitatively on the finite part of the mass shift term was already 
observed in the electric field case pl|, although then only by numerical means. 



In order to be able to do the two-loop Borel analysis, we would like to express these expansion 
coefficients as a series of T function terms. First, rewrite ( A.16| ) in the following form 
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Next, by standard combinatorics one can derive the following rearrangements, which are formally 
true for any sequence d m : 
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In these formulas, the 5 m denote the Stirling numbers of the first kind. With the help of these 
formulas, we can rewrite ( |A.17 ) as 
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where 
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Using the explicit form ( A,18| ) of dm we have found, using MATHEMATICA, the following 
remarkably complex identity (valid for m > 3): 
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Here dk are defined in ( |A.18 ). This identity has the consequence that the expansion coefficients 
in the decomposition ( [A, 21 ) are finally very simple: 
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In other words, we have shown that, up to exponentially small corrections, as n — > oo, 
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This is the expansion needed to perform the Borel analysis of the two-loop weak field expansion 
in Section g : see ( pUp and ( &Wj ). 
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